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ABSTRACT 

Coronal shocks are an important structure but without direct observations 
in solar and space physics. The strength of shocks plays a key role in shock- 
related phenomena, such as radio bursts, SEP generation and so on. This paper 
will present an improved method of calculating Alfven speed and shock strength 
near the Sun. In the method, observations as many as possible rather than one- 
dimensional global models are used. Two events, a relatively slow CME on 2001 
September 15 and a very fast CME on 2000 June 15, are selected to illustrate 
the calculation process. The calculation results suggest that the slow CME drove 
a strong shock with Mach number of 3.43~4.18 while the fast CME drove a 
relatively weak shock with Mach number of 1.90~3.21. This is consistent with 
the radio observations that a stronger and longer decameter-hectometric (DH) 
type II radio burst is found during the first event and a short DH type II radio 
burst during the second event. Particularly, the calculation results explain the 
observational fact that the slow CME produced a major solar energetic particle 
(SEP) event while the fast CME did not. Through the comparison between 
the two events, the importance of shock strength in predicting SEP events is 
addressed. 



Subject headings: Sun: coronal mass ejections (CMEs) - shock waves - particle 
acceleration 
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Introduction 



Shock waves in the corona and interplanetary (IP) medium are usually generated when 
the speeds of coronal mass ejections (CMEs) exceed the local magnetosonic wave speed. 
Shocks particularly interest researchers since they are not only associated with violent erup- 
tions from the Sun, but also may have strong capability of accelerating particles and therefore 
lead to many notable consequences which possibly influence the geo-space significantly. For 
example, gradual solar energetic particle (SEP) events, one of the major space weather effect, 
are formed as the flux of energetic particles near the Earth exceeding a threshold. Most such 
energeti c particles are thought to be generated at CME/inte rplanetary CME (ICME)-driven 
shocks (ICane et al.lll98ll ; lReameslll999l ; ICane fc Larid 120061 ). Type II radio bursts, another 
significant phenomenon in the corona and IP space, are excited by energetic electrons which 



also generally accelerated by shocks (e.g.. iNelson &: Melrose! Il985l ). Acceler ating particles is 



a dis tinct ability of shocks that could be rated by the shock strength (e.g.. Uones fc Ellison 
19911 ). Except those IP shocks arriving at 1 AU and recorded by in-situ spacecraft, many 
shocks of interest, e.g., the coronal shocks within about 5 solar radii, are far away from 
the observers and only have remote measured data. Thus how to accurately deduce shock 
strength is an inescapable step to understand the physical processes involved in shock accel- 
eration of particles and other related phenomena. 

Up to now, all coronal shocks are identified based on remote sense data. They could 
be diagnosed, for exam ple, in white-light ima ges of corona from Large Angle Spectrometric 



Coronagraph (LASCO, Brucck ner et al 



19951 ) on board Solar and Heliospheric Observatory 



(SOHO) as done by lVourlidas et al.l (120031 ). a nd in UV spectra observ ations from SOHO Ul- 



traviolet Coronagraph Spectrometer (UVCS. iKohl et al 



Mancuso et al. 2002; Raouafi et al. 2004; Ciaravella et al 



995) 



Raymond et all 12000 



20051 . 120061 ). So far, tens coronal 



shocks have been identified in UVCS spectra. However, the best signature/indicator of coro- 
nal and IP shocks is still the decameter-hectometric (DH) to kilometric type II radio bursts 
(call type lis for short) which have been known for decades. Type lis provide the informa- 
tion of plasma density at shock s, and seem to be unavoidable observations to estimation of 
shock speed and strengt h (e.g., 
Warmuth fc Mannll2005h . 



Mancuso et al. 2002; Mann et al. 2003; Vrsnak et al. 2004 



The strength of a shock could be described by the fast-mode magnetosonic Mach number 
M ms = (V s h — V sw ) /V ms , where V s h is the speed of a shock, V sw is the background solar wind 

speed, V ms = ^ v °- +c * + ^ a ~ c <*\ i s the fast-mode magnetosonic speed by assuming that the 
wave vector is along the magnetic field lines, V a is the Alfven speed , and C s is the sound 



speed. The Alfven speed is defined as V a 



B 



where B is the strength of the background 



magnetic field and p is the background plasma density. Thus the values of V s h, V sw , V a 
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and C s are required to calculate shock strength. V s h can be estimated according to the 
propagation of CMEs in LASCO field of view (FOV), which will be addressed in the next 
section. V^,, is generally about 150 km/s at 5 R s , the solar radius, and even smaller below 
that height (jSheeley et al.lll997l ). Therefore, in most cases, V sw could be ignored for a shock 
within 5 R s because the speeds of shock-associated CMEs may reach several hundreds or 
even thousands kilometers per second at 2 R s . V a involves magnetic field and plasma density, 
and C s is related to temperature. As to temperature, an isothermal atmosphere is a good 
approximation because the coronal temperature is typically at the order of one million Kelvin 
throughout the IP space. Finally, two parameters, magnetic field and plasma density, are 
left unresolved. 

Due to the lack of direct observations of magnetic field and plasma density, (quasi-)one 
dimensional global models of magnetic field and density are widely used by researchers. For 
magnetic field, the simplest 1-D model of magnetic field is described as B(r) = B s (^r) 2 , based 
on the magnetic field conservation B(r)r 2 = const, in which B s = 2.2 G is the magnetic field 
strength at solar surface R = R s . To integrate some features of active regions, the model is 
so metimes mod ified by be ing superimposed a magnetic dipole for active regions as adopted 
bv lMann etaD (e.g. l2003j)- For plasm a density, there a re also some 1- D global models, such 
as the Sait o model (ISaito et a 
them (e.g.. iGopalswamy et al. 



19771) . Newkirk model (jNewkirklll96ll ) and combinations of 



20011 ) . Combined with type II dynamic spectra, these density 



models could be used to deduce the height and speed of shocks. The fatal flaw of these models 
is also obvious. These models are highly ideal, and can resolve the variation neither in the 
plane perpendicular to the radial direction nor in temporal, which i s probably signifi cant 
as seen in section 3 for the variation of magnetic field strength and ICho et al.l (120071 ) for 
density. 

In this paper, we will present an improved method of calculating coronal Alfven speed 
and shock Mach number above 2 R s through investigating two CME events. In the method, a 
more realistic model is used to extrapolate the coronal magnetic field based on photospheric 
measurements, DH type II observations are analyzed to obtain the density at shocks, and 
CME snapshots recorded by coronagraphs are adopted to estimate the height information. 
The two events selected to illustrate the calculation process are a fast CME without a major 
SEP event and a relatively slow CME causing a major SEP event, respectively. In this 
paper, a major SEP event is defined following NOAA's criterion that the peak intensity 
of proton flux with energy >10 MeV observed at the near Earth space exceeds 10 pfu 
(particles / (cm 2 — s — ster)), which is also widely used by many other researchers. As 
will be seen, the calculation results are well consistent with observations. Particularly, the 
two studied events are right to show the importance of shock strength in predicting SEP 
events. The two events are counter examples for the picture that faster CMEs produce 
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larger SEP events. Our derived shock Mach numbers for the two events just explain the 
unusual phenomena. The next section will introduce the calculation process step by step. 
In section 3, the two events are presented to further illustrate the method. A comparison of 
the consequences between the two events with respect of SEP intensity is given in section 4. 
Finally, in section 5, we summarize the paper and make a brief discussion. 



2. Our method 



Step 1: Calculation of the local plasma density by using the dynamic spectra 
of DH type lis. Solar radio bursts of type II are indicative of s hocks propagation in 



the corona and heliosph e re accompanied by the electron acceleration (IMalitson et al.lll973 



Wild fc McCreadylll950l ; iHolman fc Pesseslll983l ). They are thought to be produced in the 
upstream plasma rather th an in the denser downstream plasma and t heir frequency reflects 
the local electron density (IDulk et al.l Il97ll ; IHolman fc Pessed Il983l ). Therefore type lis 
can be used to deduce the density of the ambient plasma that shocks are passing through. 
Type lis usually have fundamental (F) and harmonic (H) components. The fundamental 
component is emitted at the place where the wave frequency is equal to the local electron 
plasma frequency. From the frequency of the type lis fundamental (F) component, we 
can calculate the local electron number density through N e = ( iosxid 3 ) 2 ( cm ~ 3 ) • Since the 
presence of multiple species, the following relationship be tween the elec tron number density 
and the plasma density in the solar corona suggested by iPriestl (119821 ) is used to calculate 
the plasma density 



p = 1.92fim p N e = 1. 92 fj,m p 



■ fpe(Hz) 

8.98 x 10 3 



) 2 (cm- 3 ) 



where m p is proton mass with the mean molecular weight /i = 0.6 (jPriestlll982l ). 



In this paper, we study the ambient conditions and the coronal shocks only in the DH 
frequency range because (1) the DH frequency range (from about 1 to 14 MHz) corresponds 
to the heliocentric distance about 2-10 R s , where CME-driven shocks are t hought to be most 
efficient for energetic particles generation (IKahlerl 1 1 994j ; ICliver et al.ll2004l ). (2) the DH type 



lis are thought to be produced at the noses of s 



rocks, not like the metric type lis whose 



source regions are not clear (e.g., 


Gary et a 




1984: 


Maia et al. 


2000; 


Gopalswamv et a 


1997; 


Cho et al. 


2007; 


Robinson & Stewart 


1985 




Klein et al. 


1999 


Gopalswamv & Kaiser 


2002). 


The space-borne instrument, Wind/ Waves 


(Boueeret et al. 


L995h. just monitors the radio 



emissions in this frequency range, and could be used to deduce the local plasma density we 
wanted. 
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Step 2: Calculation of the heliocentric distance where DH type II radio bursts 
emitted and the speed of shocks at that time by using SOHO/LASCO images. 

SOHO/LASCO C2 and C3 cameras give the time sequence of the CME propagation from 
1.5 to 30 solar radii. It can be used to measure the height-time (HT) curve of a CME, 
and read the project ed heliocentric distance a nd s peed of the CME leading front at any 
given time. Following iGopalswamy et al.l (120051 ) and lCho et al.l (120051 ). we assume that: (1) 
DH type II radio bursts are generated at the noses of CME-driven shocks, (2) the shock 
standoff distance is relatively small near the Sun, and can be therefore ignored. For metric 
type li s, many 2-D radio images show that metric type II sources are behind the CME 
fronts (IRobinson fc Stewartl Il985l ; iKlein et al.l Il999l ). A possible explanation is that the 
shock s resulting in metric type lis are located at the flank of CMEs due to their expansion 

20021 ). However, type lis with longer wavelength are thought 



Gopalswamy fc Kaiser 



(e.g., 

to originate ahead of the CME fronts because CMEs propagate forward and shocks will be 
driven ahead of themselves. Thus the first assumption, which is also mentioned in the step 
1, is reasonable. The second assumption is also acceptable since th ere are some evidence 



that the CME bright leading edges are probably shock fronts (e.g., IVourlidas et al.l 12003 



Ciaravella et al.ll2006l ). Under these assumptions, CME leading fronts may act as a proxy of 



shocks and therefore be related to DH type lis. Besides, we assume that CMEs propagate 
radially, and correct every measured height, H measure , in 2-D LASCO images to heliocentric 
distance, H he i iocentric distance, by the equation 



H, 



measure 



heliocentric distance 



■sj (cos 6 sin <p) 2 + sin 2 6 



(2) 



where 9 is the latitude and $ is the longitude of the source region of a C ME, which is deter- 



mine d by examining SOHO/EIT (Extreme Ultraviolet Imaging Telescope, iDelaboudiniere et al. 



19951 ) movies. Then we fit the HT scatter plot with a second-order polynomial to get the real 



CME trajectory. To assess how good the approximation is, we extrapolate the fitted curve 
down to the solar surface, H heliocentric distance = 1R S , where a CME originates, and compare 
the extrapolated CME onset time Tql with the observed onset time Tqe defined as the first 
appearance of the CME-related eruptive activity in EIT images. A small difference between 
Toe and Tol indicates a good approximation since the cadence of EIT data is merely about 
12 minutes. 



Step 3: Calculation of the background magnetic field strength at the shock by 
the CSSS model. Up to now, there are no observations of 3-D magnetic field in the 
corona. Most information of coronal magnetic field comes from extra polation with the aid 
of models, such as the pot ential field source surface model (PFSS 



Schatten et al. 



Altschuler fc Newkirklll969l ). the potential field-current sheet model (PFCS. ISchatten 



1969 



19711), 
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the horizontal current-source surface model ( 
horizontal current- current sh eet model (HCCS, 



HCSS. IZhao fc Hoeksemal Il992 . |l994bj), the 



source surface model (CSSS. IZhao fc Hoeksemalll995l ; IZhao et al 



Zhao fc Hoeksema Il994af) . the current sheet 



20olx and so on. Not like 



the 1-D global models mentioned in the introduction, these models are more robust and 
complicated, in which measured photospheric longitudinal magnetic field is input as the 
bottom boundary. The extrapolated ma gnetic fields are usu ally well match the observed 
large-scale features in other wavelengths (jSchrijver et al.ll2005l ). 



In this paper, the mature CSSS model is used to extrapolate the coronal magnetic 
filed. The bottom boundary is adopted from the WSO (Wilcox Solar Observatory) synop- 
tic charts^, which are assembled from individual magnetograms observed over a Carrington 
rotation. Previous works by Zhao and his colleagues have shown that the extrapolated mag- 
netic field can appro ximately reproduce the observed variation of magnetic field in latitude 
and values at 1 AU ( Zhao fc Hoeksema 1995 ; Zhao et aD I2OO2I ) . In our study, the strength 
of the magnetic field at shocks is an average of the extrapolated magnetic field over a circle 
region with a radius of 10° at the height, where the DH type II emitted, above the CME 
source region. One may notice that the extrapolated magnetic field is a kind of average over 
a Carrington rotation and shall not exactly reflect the state at the time of interest. However, 
since what we studied is related to a large scale phenomenon and what we wanted is non- 
CME-disturbed magnetic field, we think the averaged state will not significantly distort the 
reality, nevertheless we would like to take an error of 20 percent of magnetic filed strength 
into account to make our calculation results more reliable. 



Step 4: Calculation of the Alfven speed and fast-mode magnetosonic Mach num- 
ber. From Step 1-3, the background plasmas density (p), the height and speed of the 
shock (Vgh), and the magnetic field strength (B) are deduced independently. As have been 
discussed in the third paragraph of Introduction, we will ignore the background solar wind 
speed and assume an isothermal corona with T = 1.4 MK that is corresponding to a constant 
C s of 180 km/s. Therefore, the Alfven speed and the fast-mode magnetosonic Mach number 
indicating the shock strength can be calculated. 
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Fig. 1. — The EIT and LASCO running-difference images of 2001 September 15 event. The 
dashed line indicates the direction along which the heights of the CME front were measured. 

3. Two Cases 

3.1. 2001 September 15 event 

Figure [1] shows the running difference images from SOHO/EIT and LASCO for this 
CME. The CME was first observed by LASCO/C2 at 11:54 UT on September 15, 2001. It 
looks like a bright bubble traveling along the southwest direction in the C2 FOV as shown in 
the right panel of Figure [TJ The EIT image (left panel) shows that the source region of this 
CME is at about S27W48, above which a large EUV dimming and opening of loops could be 
seen. EIT images and the soft X-ray profile from GOES satellite (the dotted line in Figure [2]) 
indicate that the CME was associated with an Ml. 5 flare. The significant eruptive signatures 
in the source region first appeared in the 11:12 UT (Toe) EIT image. By combining the 
images taken by C2 and C3 cameras, we measure the HT profile of the leading front of this 
CME until it faded into background. The direction of the measurement is chosen artificially, 
but mainly along the center of the major part of the CME as indicated by the dashed line 
in Figure [TJ The measured heights are further corrected by Equation [2J and denoted by 



http : / / wso .Stanford . edu/synopticl. html 
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cross (from C2 measurements) and asterisk (from C3 measurements) symbols in Figure [2j 
Through a second-order polynomial fitting, the heliocentric distance of the CME front as a 
function of time is obtained. The acceleration is estimated as —12.72 m/s 2 , and the CME 
speed decreased from 772 km/s at 2 R s to 674 km/s at 10 R s . From the fitting curve in 
Figure El the extrapolated onset time, Tol, of the CME is about 11:09 UT. The difference 
between Tol and Toe is only 3 minutes indicating that this fitted HT profile is a good 
approximation to the reality. 

This CME was accompanied with a clear DH type II radio burst observed by the 
Wind/Waves as exhibited in Figure [3j Since there were no other faster CMEs occurring 
near the CME of interest, the association between the DH type II and the CME is definite. 
The beginning and end times of this DH type II event have been marked by TJ, and T e in 
Figure [2j The main emission from 11:51 UT (T&) to 12:05 UT (T e ) is the fundamental, and 
the harmonic just appear at the very top of the detectable frequency range near the end of 
the event (M. L. Kaiser 2006, private communication). Since the fundamental frequency is 
diffused in a range, two white dashed lines are used to denote the upper and lower limits 
of it. At the beginning of the event, the frequency covers the range from 11.06 to 13.07 
MHz approximately, which are corresponding to the plasma density range from 2.91 xl0~ 18 
to 4.06 xl0~ 18 g/cm 3 . Both the values are used in the succedent calculations of V a , V ms , 
and M ms for showing the upper and lower limits of them. Correspondingly, the heliocentric 
distance and speed of the CME are independently estimated as 3.73 R s and 752 km/s at TJ,, 
and 4.62 R s and 741 km/s at T e . 

Following the instruction in Section 2, the coronal magnetic field extrapolation is per- 
formed by using CSSS model based on the WSO synoptic chart of Carrington rotation 1980 
corresponding to the period from 13:38 UT 2001 August 23 to 19:51 UT 2001 September 
19. Figure H] shows the strength distribution of the extrapolated magnetic field at the height 
of 3.73 R s , where the DH type II just started. It is found that the magnetic field strength 
varies in a large range from 0.03 to 0.12 G, which indicates that the Alfven speed at the 
same altitude could be largely different that directly influences the formation/strength of 
shocks. The enclosed region right above the CME source region is treated as an effective 
region, where the DH type II emitted. The mean strength of the magnetic field at this height 
in this region is 0.11 G. 

Combine the magnetic field strength with the plasma density, the Alfven speed and fast- 
mode magnetosonic speed are calculated as 124~219 km/s and 180~219 km/s, respectively, 
at the beginning of the DH type II. Therefore the shock strength, indicated by the fast- 
mode magnetosonic Mach number, of the CME is 3.43~4.18. The ranges of these values 
have considered the 20% error of the extrapolated magnetic field. Usually a shock with 
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Alfven Mach number larger than 2.7 is a super-critical shock, like the bow shock ahead of 
the magnetosphere of the Earth, while many interplanetary shocks are sub-critical ones. The 
coronal shock driven by the CME on 2001 September 15 having a fast-mode magnetosonic 
Mach number larger than 3.43 means that, although the shock driver was not so fast, the 
shock was indeed strong and had capability to cause a strong DH type II, which is factually 
observed by the Wind/Waves. Moreover, this strong shock also caused a major SEP event, 
in which the maximum proton flux intensity with energy > 10 MeV observed by GOES was 
11.6 pfu. 



3.2. 2000 June 15 event 

The CME and its source region are presented in the running-difference LASCO and 
EIT images of Figure [5j This is a fast CME first appearing at 20:06 UT on 2000 June 15 in 
LASCO/C2 and traveling along the northwest direction. It originated from N19W56, where 
a sudden brightening could be found. An Ml. 9 class flare associated with this event could 
be identified according to EIT images and the soft X-ray profile shown in Figure [6j The 
first eruptive signature recorded by EIT was at 19:36 UT. Figure [6] shows the HT scatter 
of the leading front of this CME, in which all measured heights have been corrected to the 
heliocentric distance, and the fitting curve by a second-order polynomial. By extrapolating 
the curve down to the solar surface, we find that the estimated onset time of the CME is 
19:33 UT, close to the EIT onset time 19:36 UT. This indicates that the fitted result is a 
good approximation to the reality. The acceleration of this CME was —64.54 m/s 2 , and the 
speed was 1451 km/s at 2 R s and 1176 km/s at 10 R s . 

Figure [7] displays the Wind/Waves radio spectrum for this event. A DH type II radio 
burst from about 19:52 to 19:56 UT is found. This radio burst is in the Wind/Waves Type 
II list B This DH type II radio burst is shorter but has obvious fundamental and harmonic 
components. The fundamental began at about 19:52 UT around 2.4 MHz, and the harmonic 
appeared at about 5 MHz. We could expect that the related CME-driven shock should be 
weaker than the shock driven by the first CME since the DH type II was weaker and shorter. 
The two dashed lines denote the upper and lower limits of the frequency of the fundamental 
component of this burst. At the beginning T& =19:52 UT, the frequency of fundamental 
component extended in the range from 2.24 to 2.52 MHz approximately that corresponds to 
the plasma density of 1.19~1.51 xl0~ 19 g/cm 3 according to Equation [TJ From the fitting 
curve in Figure El the heliocentric distance and speed of the CME are 3.32 R s and 1409 



2 http:/ /lep694. gsfc.nasa.gov/waves/waves. html 
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km/s at the beginning time, of the DH type II, and 3.80 R s and 1393 km/s at the end 
time, T e . It was much faster than the 2001 September 15 CME. 

The WSO synoptic chart of Carrington rotation 1963, which covers the period from 
03:06 UT 2000 May 17 to 08:17 UT 2000 Jun 13, is used to extrapolate the coronal magnetic 
filed. We use the previous synoptic chart rather than the one right covering the CME onset 
time, because the time of the CME source region passing across the central meridian is in 
Carrington rotation 1963. Figure [S] exhibits strength distribution of the magnetic field at 
H = 3.32R S , where the DH type II began to emit. The magnetic filed strength varies from 
0.033 to 0.144 G, a larger range than that in the first event. The effective region above the 
CME source region is enclosed by the black circle in Figure [H The mean magnetic field 
strength at this height in the effective region is 0.076G. 

The calculation gives the values of Alfven speed and fast-mode magnetosonic speed for 
this event 439~741km/s and 439~741km/s, respectively, at the beginning of the DH type 
II. Therefore the shock strength, M ms , of the CME is 1.90~3.21. The smaller Mach number 
indicates that the shock was relatively weak, as expected before. The lower limit of M ms , 
1.90, indicates that the shock may be as weak as that just exceed s the critical cond ition 



(Mach number is 1.6) of a shock generating suprathermal electrons (jMann et al.lll995l ). 



4. Comparison of the two cases with respect to the SEP flux 



Fast CMEs are usuall y thought to be a good pr oducer of SEP events since they usually 
drive a shock ahead (e.g., iKahlerl 1200 ll : lLariol 120051 ) . It is realized that the SEP events are 
not always fully observed near the Earth even if a CME does generate them. This is because 
these energetic particles generally move along the spiral interplanetary magnetic field lines, 
and only the particles on the field lines connecting to the Earth could be recorded. Thus 
those fast CMEs originating from western hemisph ere are thought to be mor e likelihood to 
cause an SEP event near the Earth. As shown by iGopalswamy et al.l (12004] ) . 51 out of 60 
SEP events are caused by CMEs originating from the western hemisphere. 

However, counter examples, that fast CMEs do not cause an SEP event while slow CMEs 
do, do exist. We ever studied 56 fast (projected speed > 1000 km/s) CMEs originating 
from the western hemisphere during 1997-2003, and found 12 CMEs did not cause an y 
enhancement of ener getic proton flux at the proton energy > 10 MeV Jshen et alJliooS ). 



Kahler fc Vourlidasl (120051) p resen ted 16 SEP-poor fast CMEs with speed large than 1060 
km/s. iReinard &: Andrews! (120061 ) also found 14 CMEs with speed > 1000 km/s did not 
cause an SEP event. A possible explanation is that the CME speed can not always reflect 
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the real shock strength. The comparison of flux profiles of protons with energy larger than 
10 MeV between the two events analyzed here just illustrates this issue. Of course, there 
is an alternative explanation for those slow events associated with SEPs, that the energetic 
particles do not come from shocks, but from something else like flares. Such phenomenon 
will not be addressed in this paper. 

The two events (numbered 1 and 2, respectively) have been summarized in Table [1] 
including the DH type II, CME location and speed, flare intensity, plasma density, magnetic 
field strength, Alfven speed , fast magnetosonic speed, and fast magnetosonic Mach number 
as well as the SEP intensity. Both the CMEs originated from the western hemisphere in 
the almost same longitude, and were associated with the similar flares which had the similar 
intensity (class of Ml. 5 and Ml. 9, respectively), so they are comparable. The 2001 September 
15 CME (number 1) is not a very fast CME. Its speed is 752 km/s, about half of that of the 
2000 June 15 CME (number 2). Based only on the CME speeds, we may expect that CME 
2 should produce a larger SEP event. However, the fact is that CME 2 did not produce a 
major SEP event, but CME 1 caused a major SEP event with the peak flux reaching 11.6 
pfu (as listed in the last column of Tabled]). 

The flux profiles of the protons with energy larger than 10 MeV recorded by GOES 
satellites have been presented as the dashed lines in Figure [2] and [61 The GOES instrument 
has a high background, but it will not affect our results, because GOES proton flux profiles 
are mainly used to determine whether or not there was an SEP event. Assuming zero 
scattering of protons propagating in the interplanetary space, the proton flight-time of about 
one hour from the Sun to 1 AU has been deducted from the proton flux profile. For CME 1, 
the proton flux increased soon after the CME onset and the maximum value of this proton 
flux is 11.6 pfu. For CME 2, the proton flux profile fluctuates (comparable to the background 
noise) at the bottom of Figure [6] by choosing the scale as the same as that in Figure [2J 
that suggests that there was no significant enhancement of the energetic proton flux. The 
observational facts are consistent with our calculated shock strengths, which indicate that 
CME 1 drove a stronger shock though it was relatively slow while CME 2 drove a relatively 
weak shock though it was very fast. 



Table 1: Com 

No. Date 


Darison of the 

DH Typo lis 
T b F 
UT MHz 


: 2001 Septemt 

CMEs 
Loc H V 

R s km/s 


>er 11 

Flares 


') event and 2000 June 15 event 

Derived Quantities 
P B V a Vms M ms 
10~ 18 g/cm 3 G km/s km/s 


SEP 
pfu 


1 2001/09/15 

2 2000/06/15 


11:51 11.06~13.07 
21:02 1.52~1.80 


S27W48 3.73 752 
N19W56 3.32 1409 


Ml. 5 
Ml. 9 


2.91~4.06 0.088~0.132 124~219 180~219 3.43~4.18 
0.119~0.151 0.061~0.091 439~741 439~741 1.90~3.21 


11.6 
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5. Conclusion and Discussion 

In this study, we present an improved method to calculate the Alfven speed and shock 
strength near the Sun (> 2R S ), in which as many observations as possible are used. We use 
DH type II radio bursts instead of 1-D global density model to deduce the ambient plasma 
density, use LASCO and EIT images to deduce the heliocentric distance and shock speed 
at the time of DH type lis emitting, and use WSO synoptic charts and CSSS model to 
extrapolate the coronal magnetic field. The plasma density, height and speed of shocks and 
magnetic field strength are all deduced independently. Two events, a relatively slow CME 
on 2001 September 15 and a very fast CME on 2000 June 15, are chosen to describe the 
calculation process. The calculation results suggest that the first event had a stronger shock 
with Mach number larger than 3.43 and up to 4.18, while the second event had a relatively 
weak shock with Mach number less than 3.21 and as low as 1.90. Here we would like to 
emphasize the lower limit of the Mach number of the second shock, which i s just above 



the cr itical point of a shock generating superthermal electrons as suggested by iMann et al. 



(119951 ). These results are well consistent with the radio observations that a stronger and 
longer DH type II radio burst is found during the first event and a shorter DH type II during 
the second event. 

The two events are further compared in respect of the capability of SEP generation. 
It is found that the first CME produced a major SEP event though it is relatively slow, 
while the second CME did not cause a significant enhancement of energetic protons though 
it is very fast. The fact also agrees with our calculations that the first CME drove a strong 
shock while the second CME drove a weak shock. The consistency in the comparisons of our 
calculations with both radio emissions and SEP flux implies that the calculation method we 
proposed here is reliable. According to the LASCO CME catalog^], the energy of the first 
CME was ~ 6.0 x 10 30 erg and that of the second CME was ~ 3.1 x 10 31 erg. This means 
that neither speed nor kinematic energy could correctly reflect the real strength of potential 
shocks. The importance of the accurately estimation of shock strength to the prediction of 
SEP events is therefore emphasized. 

To see how large different between our calculations and those by ideal models, Table 
[2] is compiled. For each event, the first row shows our calculated parameters that have 
been s hown in Table [H and the second row lists the plasma density by one fold Newkirk 



model (jNewkirklll96ll ). the magnetic field strength by B{r) = B s {^-) 2 with B s = 2.2 G, and 
the corresponding fast-mode magnetosonic speed and Mach number. The results by ideal 
models suggest that the first CME drove a relatively weak shock and the second CME drove 



3 http: / / cdaw.gsfc.nasa.gov/CMEJist / 
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a stronger shock. This is opposite to our conclusion, and obviously contrary to the radio 
and SEP observations shown before. 

Why the slower CME drove a stronger shock can be answered according to Table [TJ 
The ambient plasma ahead of CME 1 got based on the DH type II is much denser than that 
ahead of CME 2 while the heights of them are almost the same. In ideal density models, 
the plasma density should decrease with increasing height. However, although the height in 
event 1 is 3.73 R s higher than 3.32 R s in event 2, the plasma density in event 1 is about 27 
times of that in event 2. This is the main reason why the fast magnetosonic speed in event 
1 is much smaller than that in event 2. In addition, one may have noticed that the variation 
of density and magnetic field strength with the time and/or location changing (even in the 
same altitude) is significant. The two studied events show that the density could vary over 
about one order, and the magnetic field strength could also vary over a few orders (Fig. H] 
and [8]). Such large variations lead to the unreliability of using CME speed as a proxy of the 
shock strength. 

SEP is an important effect of space weather. In this paper, we only regard the peak 
intensity of energetic proton flux with energy larger than 10 MeV, one of the most concerned 
issue in space weather research/application. Actually, SEP acceleration and its manifestation 
are very complicated even in the events we studied here. Although there was no significant 
enhancement of the proton flux in event 2, the enhancement of the electron flux (up to 
~ 10 5 particles/cm 2 -s-ster-MeV for energy range of 38-53 keV) observed by ACE/EPAM 
instrument was much larger than that (~ 2 x 10 4 particles/cm 2 -s-ster-MeV) in event 1, and 
the enhancements of heavy ions in event 2 were also larger than those in event 1. To fully 
address these phenomena needs the knowledge of the sources of these energetic particles 
(they may come from different places, shocks and/or flares), the priorities of species to be 
accelerated under different acceleration mechanisms, and so on, which is beyond the scope 
of this paper. 

Even so, how to predict the intensities of SEP events is still an interesting enough topic. 
Although it is well known that shock acceleration is the main mechanism of SEP generation, 
the correlation between shock strength and SEP intensity is still not well established. There 

Table 2: Comparison of our results with ideal models 



No. 


Date H 


V 


P 


B 


V m s 


M ms 




R s 


km/s 


10" 18 g/cm 3 


G 


km/s 


km/s 


1 


2001/09/15 3.73 


752 


2.91-4.06 


0.088-0.132 


180-219 


3.43-4.18 




Ideal Models 




1.17 


0.16 


413 


1.82 


2 


2000/06/15 3.32 


1409 


0.119-0.151 


0.061-0.091 


439-741 


1.90-3.21 




Ideal Models 




1.62 


0.20 


442 


3.19 
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is a general re lationship bet ween CME speed and SEP intensity; however the scattering 

A major reason is that the CME speed can not reflect the 



is much large (IReames 



2000). 



real strength of the CME-driven shock. Besides, based on our calculation, the presence of 
shocks seems not to be the sufficient condition of SEP generation, there may be a threshold 
in strength, only above which the shock is able to accelerate energetic protons efficiently. 
Thus, extending our present work to a statistical study that we are doing now should be of 
benefit to our understanding of the proton acceleration at shock and space weather research. 
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Sep. 15, 2001 




10 12 14 16 18 

Time (UT) 

Fig. 2. — The heliocentric height-time plot of 2001 September 15 CME. The cross and 
asterisk symbols denote the C2 and C3 measurements, respectively. The solid line is the 
second-order polynomial fitting line. T OL denotes the extrapolated CME onset time. T OE is 
the observed CME onset time in EIT images. T b denotes the begin time of the associated 
DH type II, and H b and V b denote the height and speed of the CME front at that time. T e 
denotes the end time of the DH type II, and H e and V e are the corresponding height and 
speed of the CME front. The dashed profile shows the flux of energetic protons with energy 
larger than 10 MeV observed by GOES. A one-hour flight time of protons from the Sun 
to 1 AU has been deducted from the proton flux profile by assuming zero scattering in the 
interplanetary space. The dotted line denotes the GOES soft X-ray flux at the wavelength 
of 1.0 - 8.0 A. 
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Wind Waves RAD2 receiver: 2001/9/15 
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Fig. 3. — The Wind/Waves observations during 11:00 - 13:00 UT on September 15, 2001. 
The dashed lines denote the frequency range of the fundamental component of a DH type 
II. 



WSO--1980 at altitude of 3.73 R s 




Carrington Longitude 

Fig. 4. — The distribution of the extrapolated magnetic field strength at the height of 3.73 
R s , where the DH type II began to emit. 
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Fig. 5.- 



The EIT and LASCO running-difference images of 2000 June 15 event. 
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Wind Waves RAD2 receiver: 2000/6/1 5 
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Fig. 7.— The Wind/Waves observations during 19:00 - 21:00 UT on June 15, 2000. 




Fig. 8. — The extrapolation magnetic field strength at the height of 3.32 R S} where the DH 
type II began to emit. 



